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Abstract 
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1 Introduction 


The reasons for interests in the specialized studies on P5-induced reactions are modestly well-known 
and are, by now, somewhat commonplace for which we are not going to elaborate on them. In the 
recent past the WA98 Collaboration P presented a detailed study of neutral pion transverse mass 
spectra in the range ^.hGeVj(? < mx — ttiq < \S^GeVjc? and 2.3 <y< 3.0 for collisions of PbPb and 
PbNb at 158A GeV for different centralities. Besides, Aggarwal et al (WA98 Collaboration) PP also 
pointed out how the data on the ratios of invariant multiplicities in nucleus-nucleus reactions to those 
in proton-proton reactions diverge from values predicted by nearly all the important standard versions 
of multiple production of particles in nucleus-nucleus collisions. Against this perspective, we would 
try to understand here the reported features of neutral pion production in two lead(A = 208)-induced 
reactions at 158A GeV and also try to explain the problematical ratios in the light of a combinational 
approach The present study is conducted obviously with the clear motivation of testing this newly 
projected approach for understanding the characteristics of particle production in heavy ion collisions. 

We wonld concentrate here, in the main, on the following few aspects: (i) the nature of pionic 
inclusive cross-section, (ii) the property of average transverse momentum, (iii) the qualitative character 
of the multiplicity-behaviour and of their ratios, (iv) the nature of centrality-dependence of the ratios 
in the various p'p-intervals. Our results on these observables and the available data would also be 
compared with some of the most prominent models for multiparticle production in heavy ion collisions. 
The paper would be organized as follows: The section 2 offers the outline of the model and the sketch 
of the physical ideas which prompted us to proceed in the stated direction. In section 3 we present 
the calculational results and their graphical descriptions. The last section is, as usual, for our final 
comments on the totality of the work presented here. 

2 The Model: A Sketch 

Following the suggestion of Faesslerj^ and the work of Peitzmann^ and also of Schmidt and Schukraft [3] , 
we propose here a generalized empirical relationship between the inclusive cross-section for pion pro¬ 
duction in nucleon(N)-nucleon(N) collision and that for nucleus(A)-nucleus(B) collision as given below: 

{AB ttX) ~ {PP ttX) , (1) 

where 4>{y, pr) could be expressed in the factorization form, (p{y, pt) = f{y) g{pt)', and the 
product, AB on the right hand side of the above equation is the product of mass numbers of the two 
nuclei participating in the collisions at high energies, of which one will be the projectile and the other 
one the target. 

While investigating any specific nature of dependence of the two variables (y and pt), either 
one of these is assumed to remain constant. In other words, more particularly, if and when the pt- 
dependence is studied by experimental groups, the rapidity factor is treated to be constant and the 
vice-versa. So, the formula for studying the nature of pr-spectra turns into 

E^ {AB ^ ttX) ~ E^ {PP ^ ttX) , (2) 

dp-^ dp-^ 

The main bulk of work, thus, converges to the making of an appropriate choice of form for g{pT)- 
And the necessary choices are to be made on the basis of certain premises and physical considerations 
which do not violate the canons of high energy particle interactions. 
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The expression for inclusive cross-section of pions in proton-proton scattering at high energies in 
Eqn.(2) could be chosen in the form suggested first by Hagedorn|H]: 

(PE ^ vrX) = Cl ( 1 + , (3) 

dp'^ Po 

where Ci is the normalization constant, and Po, n are interaction-dependent chosen phenomeno¬ 
logical parameters for which the values are to be obtained by the method of fitting. 

The final working formula for the nucleus-nucleus collisions is now being proposed here in the 
form given below: 


{AB ttX) oc (AB)^^ + - dvl) {PP nX) 

oc {AB)<^^ + - dpi) (1 + ^ 

with g{pT) = (e -|- apT — PPt)^ where this suggestion of quadratic parametrization for g{pT) is 
exclusively made by us and is called hereafter De-Bhattacharyya parametrization(DBP). In the above 
expression e, a and (3 are constants for a specific pair of projectile and target. 

Earlier experimental works^EllHl showed that g{pT) is less than unity in the pT-domain, pT < 1.5 
GeV/c. Besides, it was also observed that the parameter e, which gives the value of g{pT) at Pt = 0, 
is also less than one and this value differs from collision to collision. The other two parameters a 
and (3 essentially determine the nature of curvature of g{pT)- However, in the present context, precise 
determination of e is not possible for the following understated reasons: 

(i) To make our point let us recast the expression for (4) in the form given below: 

E^{AB ttX) ^ C 2 {ABf {AB)^^P^-dpl) (1 (5) 

dp-^ Po 

where C 2 is the normalization term which has a dependence either on the rapidity or on the rapidity 
density of the pion and which also absorbs the previous constant termjCi as well. 

Quite obviously, we have adopted here the method of fitting. Now, in Eqn.(5) one finds that 
there are two constant terms C 2 and e which are neither the coefficients nor the exponent terms of 
any function of the variable, px- And as e is a constant for a specific collision at a specific energy, 
the product of the two terms C 2 and (A.BY appears as just a new constant. And, it will just not 
be possible to obtain fit-values simultaneously for two constants of the above types by the method of 
fitting. 

(ii) From Eqn.(2) the nature of gipx) can easily be determined by calculating the ratio of the 
logarithm of the ratios of nuclear-to-PP collision and the logarithm of the product AB. Thus, one 
can measure e from the intercept of gipx) along y-axis as soon as one gets the values of for both 
AB collision and PP collision at the same c.m. energy. But, there is a problem that it will not be 
possible to get readily the data on inclusive spectra for PP collisions at all c.m.energy values. 

In order to sidetrack these difficulties and also to build up an escape-route, we have concentrated 
here almost wholly to the values of a and (3 for various collision systems and the resultant effects of 
C 2 and e have been absorbed into a single constant term Co- Hence, the hnal expression becomes 


E^(AB ^ ttX) ^ C 3 (AB)^^P^-dPl) (1 +^)-^ 
dp'^ Po 


( 6 ) 


with C 3 = C 2 {ABY- 

The exponent factor term apx — PPt obviously represents here [g^px) — c] instead of g{px) alone. 
The expression(6) given above is the physical embodiment of what we have termed to be the grand 
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combination of models(GCM) that has been utilized here. The results of PP scattering are obtained 
in the above on the basis of eqn.(3) provided by Hagedorn’s model(HM); and the route for converting 
the results of NN to NA or AB collisions is built up by the Peitzmann’s approach(PA) represented 
by expression(2). The further input is the De-Bhattacharyya parametrization for the nature of the 
exponent. Thus, the GCM is the combination of HM, PA and the DBP, all of which are used here. 

And the choice of this form of parametrization for the power of the exponent in eqn.(4) is not 
altogether a coincidence. In dealing with the EMC effect in the lepton-nucleus collisions, one of the 
authors here(SB),jn] made use of a polynomial form of A-dependence with the variable Xir(Feynman 
Scaling variable). This gives us a clue to make a similar choice for both g{pT) and f{y) variable(s) 
in each case separately. In the recent times, De-Bhattacharyya parametrization is being extensively 
applied to interpret the measured data on the various aspects (10] of the particle-nucleus and nucleus- 
nucleus interactions at high energies. In the recent past Hwa et. al. El also made use of this sort of 
relationship in a somewhat different context. The underlying physics implications of this parametriza¬ 
tion stem mainly from the expression(4) which could be identified as a clear mechanism for switch-over 
of the results obtained for nucleon-nucleon(PP) collision to those for nucleus-nucleus interactions at 
high energies in a direct and straightforward manner. The polynomial exponent of the product term 
on AB takes care of the totality of the nuclear effects. 

For the sake of clarity and confirmation, let us further emphasize a point here very categorically. 
It is to be noted that this model(GGM) containing all the Eqns.(4), (5) and (6) was described in some 
detail earlier and was made use of in analyzing extensive sets of data in the previous pubhcatioins|2lEj 
by the same authors. And in verifying the validity of this model further, the purpose here is to apply 
the same model to some other problematical aspects of data which we would dwell upon in the 
subsequent sections. 

Indeed, quite obviously, there are two phenomenological parameters in g{pT) which need to 
be physically explained and/or identified. In compliance with this condition we offer the following 
physical explanations for the occurrence of all these factors. The particle-nucleus or nucleus-nucleus 
collisions at high energies almost instantly gives rise, ex hypothesi, to an expanding blob or fireball 
with rising temperature. In real and concrete terms this stage indicates the growing participation of 
the already-expanded nuclear blob. As temperature increases at this stage, the emission of highly 
energetic secondaries (which are mostly peripheral nucleons or baryons) with increasing transverse 
momentum is perfectly possible. The coefficient a addresses this particularity of the natural event; 
and this is manifested in the enhancement of the nuclear contribution with the rise of the transverse 
momentum. Thereafter, there is a turnabout in the state of reality. After the initial fractions of 
seconds, the earlier-excited nuclear matter starts to cool down and there is a clear natural contraction 
at this stage as the system suffers a gradual fall in temperature. Finally, this leads to what one might 
call ‘freeze-out’ stage, which results in extensive hadronization, especially in production of hadrons 
with very low transverse momentum. In other words, the production of large-pr particles at this 
stage is lowered to a considerable extent. This fact is represented by the damping or attenuation 
term for the production of high-p'r particles. The factor (5 with negative values takes care of this 
state of physical reality. Thus the function denoted by g{pT) symbolizes the totality of the features 
of the expansion-contraction dynamical scenario in the after-collision stage. This interpretation is, at 
present, is only suggestive. However, let us make some further clarifications. 

The physical foundation that has here been attempted to be built up is inspired by thermody¬ 
namic pictures, whereas the quantitative calculations are based on a sort of pQCD-motivated power- 
law formula represented by eqn.(4). This seems to be somewhat paradoxical, because it would be 
hard to justify the hypothesis of local thermal equilibrium in multihadron systems produced by high 
energy collisions in terms of successive collision of the QCD-partons(like quarks and gluons) excited 
or created in the course of the overall process. Except exclusively for central heavy ion collisions, a 
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typical parton can only undergo very few interactions before the final-state hadrons ‘freeze out’, i.e. 
escape as free particles or resonances. The fact is the hadronic system, before the freeze-out starts, 
expands a great deal - both longitudinally and transversally - while these very few interactions take 
place [T^. But the number of parton interactions is just one of the several other relevant factors for the 
formation of local equilibrium. Of equal importance is the parton distribution produced early in the 
collision process. This early distribution is supposed to be a superposition of collective flow and highly 
randomized internal motions in each space cell which helps the system to achieve a situation close to 
the equilibrium leading to the appropriate values of collective variables including concerned and/or 
almost concerned quantities. The parameter a in expression(4) is a measure of the ratio of the net bi¬ 
nary collision number to the total permissible number among the constituent partons in the pre-freeze 
out expanding stage identified to be a sort of explosive ‘detonation m stage. This is approximated 
by a superposition of collective flow and thermalized internal motion, which is a function of hadronic 
temperature manifested in the behaviour of the average transverse momentum. The post freeze-out 
hadron production scenario is taken care of by the soft interaction which is proportional ^ E] to the 
number of participant nucleons, Npart, according to almost any variety of wounded nucleon model. 
The factor (5, we conjecture, offers a sort of the ratio of the actual participating nucleons to the to¬ 
tal number of maximum allowable(participating) nucleons. In fact, this sort of physical explanations 
seems to have been acceded to by some of the physics community through their approval of some of 
our previous works [THl I14j . 

3 The Calculations and the Results 

The graphical plots presented in the diagrams in Fig. ^ to Fig. |S1 describe the measured data on pion 
production modestly well in the region of the moderate values of the transverse momentum. The Fig.l 
presents the fit for data on PP collisions. As described in the previous section the parameter values 
to be inducted in calculations of PP cross sections are to be obtained from Table-1. These are based 
on Hagedorn’s model. Besides, in order to arrive at the theoretical values of invariant pion production 
cross section in PbPb and PbNb collisions at the measured energy with the help of this combinational 
approach, one has got to use the fit values of a and j3 as provided in Table-2 and Table-3 for minimum 
bias events and for the various centralities of the nuclear collisions respectively. The average number of 
binary nucleon-nucleon collisions, < Ncoii > and that of participant nucleons, < Npart > at different 
centralities are to be obtained from Table-4 which is an adaptation of a specific set of simulation 
resultspp. The statement on the nature of agreement with regard to the invariant multiplicity shown 
in Fig. 01 for pion production in PbPb collisions at 158A GeV in the most peripheral and most central 
collisions, and also in the minimum bias events remains valid. These results obtained by the GCM 
are compared in Fig. |1] with some of the very popular models in the field, while the comparisons 
with the extracted datapP as well lay in the background in all cases. Compatibility with data over 
a wide range of the p^-values is modestly satisfactory. The graphs on the nature of the average 
transverse momentum vs. the number of participating nucleons which is a measure of the centrality of 
the collisions, depict faithfully (Fig. the expected behaviours. The average transverse momentum 
values here are normally defined by: 


<PT > 


/p Jin PT dpT 

rPr dN j 

Jpmin iiprj^ UPT 


(7) 


where p™*"" and are the lower and upper limits of the transverse momentum, pp over which 

the integration has been performed. The values of average transverse momenta, < px > have been 
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calculated in three rriT—mQ intervals; and they are 0.5 GeV/(? < rriT—mQ < 2.0 GeV/c^, 1.0 GeV/c^ < 
mT — niQ < 2.0 GeV/(? and 2.0 GeV!(? < rriT — fnQ < 3.0 GeVjc?. Here mr and tuq are the transverse 
mass and the rest mass of secondary pions respectively; and they are related with p-p by the equation 

rn}p = Pj' + rriQ. 

However, for a specific nucleus-nucleus collision, the nature of dependence of any observable on 
the center of mass(c.m.) energy per basic colliding pair of nucleons, i.e.^SNNj is to be obtained from 
the same for what it is in nucleon-nucleon(PH) reactions at high energies. 

The ratio-behaviours of the pr-spectra and the nature of another observable as would be defined 
later are studied here in the rest of the diagrams. The ratio-values are of two types. In Fig. |H1 
the ratio-values of the inclusive p^-spectra for various degree of centrality of collisions to the same 
for minimum bias events have been studied. The other set of ratios are for the multiplicities for 
various centrality-degrees of collisions. The diagrams shown in Fig. |7(a)]|7(d)| demonstrate the plots of 
the ratios of invariant multiplicity distributions of neutral pions normalized to the number of binary 
collisions versus the transverse momentum. In Fig. |7(a)| the ratio of the invariant cross section of 
neutral pions produced in least central PbPb collision(A^( 7 o;; ~ 10) normalized by the number of 
binary collisions to that of neutral pions produced in PP collisions has been plotted as a function of 
transverse momentum; and the ratio-values, in this particular case, are given by 


Ratio{ 


Pb + Pb 
P + P ^ 


d'p'^ I 



( 8 ) 


In the figure(Fig. |7(a)| ) the extracted data shown by the filled squares were obtained from Aggar- 
wal et alpP who assumed some parametrization for PP reaction. In case of GCM, in order to obtain 
the values of E^^\pp, we have first utilized the parameter values given in Table-1 for the E^^\pp 
and then normalized them by the total inelastic cross section ain = 43.6 ± 4.0 mb CS], according to 
the conventional rule of conversion arising out of the definitions. 

The diagrams in Fig. |7(b)]|7(d)| depict the ratios of invariant cross sections at different centralities 
with a normalization of the corresponding number of binary collisions. The ratios in these cases can 
be written in the form 


Ratio{ 


Ncoii = X 
Ncoii = y 


dp^ I 

NcoU 

p d^N 
dp^ I 

Ncoii 


(9) 


where X and Y denote two different number of binary collisions at two different centralities. 
The observable plotted in all the diagrams of Fig. |8(a)||8(d)| and denoted by Rbim is defined by 
the relation , 


Rbin = E-j^/Ncoii ( 10 ) 

The extracted datapP on the observable related to neutral pion yield per binary collisions are for 
different centralities in the separate pj’-intervals. 

Both the sets of figures(Fig. |7(a)||7(d)| and Fig. |8(a)|8(d)| ) provide roughly satisfactory description 
of the data. So, the present work could be presented as a successful continuation of one of our previous 
work [2]. Thus all this lends substantial degree of credence to the combinational approach built up 
and completed by us. And in comparison with even the most accepted version of the simulation-based 
standard models, i.e. the HIJING model, the present approach works quite well. 
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4 Concluding Remarks 


Based on the analyses made by the chosen approach here, we make the following particular observa¬ 
tions: 

(i) The modestly successful reproduction of the measured data on PP reaction, PbPb and PbNb 
collisions in the minimum bias events confirm that the used basic models have good degree of 
reliability. 

Besides, the efficacy of the model is exhibited in a somewhat tangible way for the various degree 
of centrality of nucleus-nucleus collisions demonstrated by Fig.|21 In fact, the Fig.EJis a summary 
statement in favour of the chosen model. 

(ii) The average px behaviour is also explained by the present approach. But for any specific nucleus- 
nucleus collision, the s-dependences of the < px > values of the various secondaries are to be 
obtained mainly from the analysis of nucleon-nucleon collision in a model-based manner. 

(iii) The ratios of the collisions with various degrees of centrality of collisions to the minimum bias 
events have been obtained with a fair degree of consistency. 

(iv) The ratios of the most central to semi-central or least central(peripheral) nuclear collisions have 
also been computed by the used approach in conformity with the extracted available datapP on 
pion production alone. 

But, the above points cover no comments on the comparison of the performances of the few chosen 
models which are under consideration here and offer no insights into the state of physical reality. 
Regarding the former, the Fig.7 illustrates somewhat convincingly that the GCM performs much 
better to accommodate data than both the PQCD and HIJING versions of the standard multiparticle 
production models separately do. And, concerning the latter, firstly, the status of anomalous nuclear 
enhancement(ANE) is still ambiguous. In so far as pion production in high energy PbPb collision is 
concerned, this effect is seen (Fig.7b) to be perceptible for the relatively peripheral collisions. And for 
these non-central collisions the GCM describes the observed effect modestly well. But the scenarios 
are different for the other sets of the medium central to most central collisions(Figs.7c and 7d), 
wherein the enhancement effect ceases to exist and the reversals of nature occur quite prominently with 
observations of the gradual and puzzling diminution of the effects related with transverse mass. Thus, 
in any of such cases, the concept of scaling with system size is not corroborated by the extracted data 
obtained by WA98 collaboration^. True, all these data-trends are captured by the GCM somewhat 
automatically. But the reasons for this spectacular success in terms of collision dynamics, i.e. the 
number of collisions or the number of participants, are still not very clear to us. Secondly, yields of 
neutral pions per binary collisions could increase with the number of collisions only upto a certain 
degree of ‘hardness> 1 GeV/c - 2 GeV/c) of the multiple collisions. For large pT-i'anges, the 
validity of the physics of parton saturation at larger transverse momentum region is quite manifest at 
values of Ncoii ~ 200(Figs.8c and 8d) from both the extracted data-sets and GCM-based results. 
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Table 1: Different parameter values for vr^ production in P + P collision at 200 GeV 


Cl 

Po(GeV/c) n 

270 ±5 

3.0 ±0.2 23 ±1 


Table 2: Necessary parameter values for neutral pion production in the minimum bias Pb + Nb and 
Pb + Pb collisions at 160A GeV. 


Gollision 

Cs 

a(GeV/c)-^ 

/3(GeV/c)-^ 

Pb + Nb 

(3.0 ±0.5) X 10^ 

0.17 ±0.02 

0.033 ± 0.002 

P6±P6 

(8.6 ±0.8) X 10® 

0.17 ±0.02 

0.032 ± 0.001 


Table 3: Various parameter values for vr^-production in Pb + Pb collisions at 160A GeV for different 
centrality-values. 


Ncoii 

Cs 

a(Gi 

eV/c)-^ 

/3(GeV/ 

c)-^ 

10 

19 ±1 

0.13 

± 

0.02 

0.042 

± 

0.002 

30 

33 ±2 

0.15 

± 

0.02 

0.035 

± 

0.003 

78 

130 ±4 

0.14 

± 

0.02 

0.038 

± 

0.003 

207 

200 ±6 

0.17 

± 

0.01 

0.036 

± 

0.003 

408 

250 ±3 

0.16 

± 

0.02 

0.030 

± 

0.003 

570 

495 ±5 

0.16 

± 

0.01 

0.032 

± 

0.004 

712 

746 ±3 

0.14 

± 

0.02 

0.035 

± 

0.003 

807 

705 ±8 

0.15 

± 

0.01 

0.040 

± 

0.003 


Table 4: Used values of < Npart > and < Ncou > for various centrality-classes pQ of PbPb collisions. 
The values in column ‘Glass’ indicate vertically downwards gradual transitions of the collisions from 
the-lowest-to-the-highest centrality. 


Glass 

PT(GeV) 

^ ^part ^ 

< Ncoll > 

1 

< 24.35 


12 ±2 

9.9 ±2.5 

2 

24.35 - 

55.45 

30 ±2 

30 ±5 

3 

55.45 - 

114.85 

63 ±2 

78 ±12 

4 

114.85 - 

- 237.35 

132 ±3 

207 ± 21 

5 

237.35 - 

- 326.05 

224 ± 1 

408 ± 41 

6 

326.05 - 

- 380.35 

290 ±2 

569 ± 57 

7 

380.35 - 

- 443.20 

346 ± 1 

712 ± 71 

8 

> 443.20 

380 ± 1 

807 ± 81 
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Figure 1: Plot of as a function of transverse momentum, px for production of secondary neutral 
pions in PP collision at E^ab = 200 GeV. The experimental data are taken from m- The solid curve 
represents the fit obtained on the basis of eqn.(3). 
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Figure 2: The inclusive spectra of secondary neutral pions produced in two induced reactions at 
Eiab = 160A GeV(minimum bias). All the experimental data are taken from pQ. The GGM-based 
results(eqn.(6)) are shown by the solid curves. 
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Figure 3: Invariant multiplicities of produced in Pb + Pb collisions of different centralities at 
Eiab = 160A GeV as a function of rriT — mo- The experimental data-points for various centrality bins 
are taken from The solid curves provide the GCM-based fits(eqn.(6)). 
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Figure 4: A comparison of the performance of different models via the ratio of the minimum bias data 
for Pb + Pb collision to model-based results. The solid squares provide the GGM-based results. Other 
model-based results are taken from 
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Figure 5: Plot of < pr >, obtained on the basis of the present model(GCM), as a function of the no. 
of participant nucleons, Npart- 
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Figure 6; Plot of ratios of invariant multiplicities of neutral pions for Pb + Pb collisions of different 
centralities to minimum bias events. The data-type points represent the experimentally measured 
values while the solid curvilinear lines show the GCM-based results. The dashed curves are for 
uncertainties arising out of the error-ranges of all the parameter values. 
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Figure 7: Ratios of of neutral pious for different centralities normalized to the number of binary 

collisions. The GCM-based results are depicted by the solid curves. The dashed curves represent the 
uncertainties arising due to the errors in Cs, a and f3. Other model-based results have been obtained 
from pp. 
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Figure 8; Invariant multiplicities of neutral pions normalized to the number of binary colli- 
sions(eqn.(10)) as a function of Ncoii for different rriT — mo intervals. The open squares represent 
the experimentally extracted values, while the fully filled ones indicate the GCM-based results. 
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